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INTRODUCTION 


It  has  been  known  for  many  years  that  nitrocellulose  (NC)  and  nitroglycerin  (NG) 
are  not  stable  and  decompose.  However,  very  little  is  known  about  the  decomposition 
piocess.  it  was  found  empirically  that  the  addition  of  certain  compounds  apparently 
slowed  the  decomposition  process  or  at  least  hindered  a  catastrophic  explosion.  This 
hindering  of  the  explosive  process  occurs  by  the  removal  of  N02  by  chemical  reaction 

with  the  additional  compounds  or  stabilizers  as  they  are  called.  It  seems  apparent  that 
diffusion  must  occur  in  these  materials  for  the  stabilizer  to  be  able  to  do  its  work,  yet 
little  or  nothing  is  known  about  diffusion  in  NC  or  NC/NG  base  propellants.  A  general 
discussion  of  diffusion  is  followed  by  a  discussion  of  an  explicit  model  for  stabilizer 
reaction  with  N02  via  diffusion  of  the  latter. 

DIFFUSION 

There  are  several  types  of  experimental  evidence  which  indicate  that  some  sort  of 
diffusion  takes  place  during  thermal  decomposition  of  propellant.  Gas  is  given  off 
during  thermal  decomposition.  Several  different  gases  have  been  detected  in  this 
process.  This  indicated  a  general  breakup  of  nitrocellulose  and  nitroglycerin  (ref  1). 
However,  it  is  not  clear  whether  these  gases  come  from  surface  layers  or  from  the  bulk, 
thus  indicating  bulk  diffusion.  Perhaps  more  direct  evidence  for  bulk  diffusion  comes 
from  the  conversion  of  the  stabilizer  into  its  nitroderivatives  by  the  addition  of  nitro 
groups  (ref  2).  It  appears  as  if  this  process  can  only  take  place  by  NOx  groups  (where  x 

is  1 ,  2  or  3)  breaking  off  from  NC  (and/or  NG)  and  diffusing  until  captured  by  the  stabi¬ 
lizer  or  its  derivatives.  Alternately,  it  is  necessary  to  assume  that  the  much  larger 
stabilizer  molecules  are  mobile  and  diffuse.  The  latter  seems  unlikely. 

Evidence  has  been  presented  for  NG  diffusion  under  some  conditions.  NG  has 
been  observed  to  collect  on  the  surface  during  heat  treatment,  thus  indicating  NG 
diffusion  from  the  bulk  to  the  surface  (ref  3).  The  exact  conditions  under  which  NG 
diffusion  takes  place  have  not  been  determined. 

The  activation  energies  for  diffusion  of  any  type  in  NC  base  propellants  have  not 
been  determined  or  identified  as  such.  Overall  or  global  activation  energies  have  been 
determined,  but  the  parts  of  these  energies  which  could  be  due  to  diffusion  have  not 
been  separated  out.  In  fact,  throughout  the  literature,  the  activation  energies  which 
have  been  measured  have  been  associated  with  bond  breaking  only  and  not  with 
diffusion.  In  simple  ionic  crystals,  such  as  alkali  halides  and  silver  halides,  very  clever 
techniques  have  been  used  to  separate  the  activation  energies  for  diffusion  from  the 
energies  for  bond  breaking.  Unfortunately  the  techniques  used  in  the  simple  ionic 
crystalline  solids  can  not  be  used  in  propellant  polymer  solids. 
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In  general,  the  global  activation  energies  which  have  been  observed  during 
propellant  thermal  decomposition  have  been  found  to  increase  with  increasing  tempera¬ 
ture  of  observation,  thus  indicating  the  possibility  of  different  processes  in  different 
temperature  ranges  (refs  1  and  4),  It  is  also  necessary  to  consider  the  possibility  that 
the  simple  gas  molecules  which  are  detected  on  thermal  decomposition,  e.g.,  NO,  N20, 

C02,  N00,  CO  and  N2  move  easily  through  the  propellant  matrix  at  least  in  some  direc¬ 
tions.  There  will  be  a  specific  activation  energy  for  diffusion  for  each  type  of  gas 
molecule  and  also  most  probably  an  activation  energy  for  bond  breaking  for  each 
molecular  species.  Clearly  more  work  must  be  done  to  establish  the  role  of  diffusion  in 
the  thermal  decomposition  of  NC  base  propellant  materials. 

Diffusion  in  polymer  systems  (such  as  NC)  is  more  complex  than  in  many  other 
types  of  solid  materials.  For  example,  the  diffusion  coefficient  may  not  be  a  constant  at 
any  given  temperature  but  may  depend  on  quantities  such  as  the  concentration  (ref  5). 
In  addition,  the  diffusion  coefficient  may  be  different  along  the  direction  of  the  polymer 
chain  tnan  perpendicular  to  it.  The  diffusion  of  many  simple  gases  has  been  measured 
in  a  variety  of  polymers.  The  diffusion  coefficients  have  been  found  to  vary  over  several 
orders  of  magnitude  for  the  same  simple  gas  molecules  in  different  polymers  (ref  6). 
Thus,  it  is  very  difficult  or  impossible  to  estimate  the  magnitude  of  the  diffusion  coeffi¬ 
cients  of  the  same  simple  gases  in  NC  systems  although  these  same  simple  gases  are 
evolved  during  thermal  decomposition. 

An  effort  was  made  to  estimate  the  maximum  number  of  thermally  activated  jumps 
an  NOx  unit  must  make  after  bond  scission  to  arrive  at  a  stabilizer  molecule  for  a  double 

base  rocket  propellant.  Many  simplifying  assumptions  were  made  including  uniform 
distribution  of  the  various  ingredients  in  the  propellant  and  also  that  an  NG  molecule 
occupies  the  same  space  as  one  NC  ring.  The  separation  of  the  stabilizer  molecules  in 
NC  ring  lengths  was  then  estimated.  The  maximum  distance  which  an  NOx  unit  must 

make  is  then  one  half  of  the  separation  distance  between  stabilizer  molecules.  The 
number  of  jumps,  n,  in  a  random  walk  process  is  given  by 


where  R  is  the  net  displacement  in  n  jumps  and  r  is  the  length  of  one  jump  (ref  7). 
Using  this  formula,  the  maximum  number  of  jumps  in  terms  of  NC  ring  lengths  which  an 
NOx  species  must  make  after  scission  from  NC  to  arrive  at  a  stabilizer  molecule  was 

fuund  to  be  about  10.  If  the  elementary  jump  distance  is  smaller  than  one  NC  ring 
length  the  number  of  jumps  will  be  larger. 
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KINETIC  MODEL 


As  an  example  of  diffusion  in  a  propellant,  consider  the  thermal  degradation  of  a 
double  base  rocket  propellant  containing  two  nitro  diphenylamine  (2NDPA)  as  the 
stabilizer.  It  is  assumed  that  the  NOx  species  produced  by  thermal  degradation  of  NC 

and  NG  diffuse  randomly  until  captured  by  the  stabilizer  or  its  derivatives.  The  latter  are 
2NDPA  +  b(NOx)  where  b  is  a  small  integer  indicating  the  number  of  NOx  captured 

previously  by  the  stabilizer  molecule.  While  it  is  not  known  whether  it  is  NO,  N02,  or 
N03,  which  diffuses,  the  net  result  of  capture  at  the  stabilizer  (or  its  derivatives)  is  the 
addition  of  a  N02  group  to  the  stabilizer  or  its  derivatives. 


For  the  purpose  of  generating  a  kinetic  model  to  describe  the  capture  of  NOx  by 
the  stabilizer  or  its  derivatives,  the  following  are  defined: 

Let  w  =  concentration  of  2NDPA 

2,4  Dinitro  DPA 

Let  y  =  concentration  of 

2,4'  Dinitro  DPA 
2,4,4'  Trinitro  DPA 

Let  z  =  concentration  of 

2. 2’, 4  Trinitro  DPA 

Let  r  =  concentration  of  2,2',  4,4’  Tetranitro  DPA 

Let  s  =  concentration  of  end  product  (simplifying  assumption  to  demonstrate 

model) 

where  all  compounds  in  the  right  column  (except  2NDPA)  are  derivatives  of  2NDPA. 
For  simplicity  it  is  assumed  that  2,4  Dinitro  DPA  and  2,4’  Dinitro  DPA  and  the  two  trinitro 
DPAs,  respectively,  have  the  same  cross  sections  for  capture  and  so  are  not  separated. 
That  this  is  approximately  the  case  can  be  seen  by  an  inspection  of  figure  1.  The 
following  reactions  are  assumed  to  occur  as  2NDPA  and  its  derivatives  capture  NOx 

groups 
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ki 

w  +  N0X  - -  y 

^2 

y  +  N0X  - ►  z 

kg 

z  +  N0X  — - — ►  r 

l<4 

r  +  NOx  - -s  (1) 

where  k, ,  k2 ,  k3 ,  and  k4  are  rate  constants  as  used  in  equations  2  through  6.  The 
differential  equations  which  describe  the  reactions  can  be  written  as 


™  =k1(NOx)w  =  k,w  (2) 

^  =k1(NOx)w-k2(NOx)y»k1w-k2y  (3) 

=k2(NOx)y-k3(NOx)z  =  k2y-k3z  (4) 

%  =k3(NOx)z-k4(NOx)r  =  k3z-k4r  (5) 

=k4(NOx)r  =  k4r  (6) 


where  the  concentration  of  diffusing  (NOx)  is  assumed  to  be  independent  of  time  to 
obtain  the  right  side  o?  equations  2  through  6.  The  boundary  conditions  are 


at  t  =  o 


w  =  wo  (7) 

y=z=r=s=o 
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The  solutions  to  equations  2  thorugh  6  with  the  boundary  conditions  are  as  follows 


w  =  w0e 


kl  Wr>  /  kl  I  k't\ 

°_  (e  e  ‘  ) 


v  -  <k2-k, ) 


k1k2W0r(ekl'-ek’') 


(k2 -ki )  1_  ( K3  -  ki)  (k3  -  k2) 


k,k2k3w0  [  (e  k,t-e k4t) 


(e  k,t-e  kdt) 


r  (k2  -  ki)  j  (k4  -  ki)  (k3  -  ki)  (k4  -  k2)  (k3  -  k2)  + 


(e  M  -e  M)  (k2  -ki) 

(k4  -  k3)  (k3  •  k2)  (k3  -  ki) 


_  _  ki  ^2  k3  k4  w0  •: _ 1 _ f  (1  -  e  k4t)  _  (1  -  e  M)  \ 

(ki  -  k2)  (ki  -  k4)  (ki  -  k3)  l  ki  ) 


(ki  -  k2)  (ki  -  k4)  (ki  -  k3)  y  k4 

L 

1  f  (l  -  e k?1)  (l-e  M) 

(k2  -  k4)  (k2  -  k3)  y  k2 


(k,  -  k2)  f  (1  -  eM)  (I-  e kj,)>  ' 

(k3  -  k4)  (k2  -  k3)  (ki  -  k3)  y  ^3  J 


In  addition  to  satisfying  the  boundary  conditions  these  solutions  have  the  following 
properties 

1 .  at  t=0  w  and  y  have  finite  slopes  (1 3) 


2.  at  t  =  0  z,  r  and  s  have  zero  slopes 
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,z  and  r  have  maxima  as  a  function  of  time 


(15) 


4.  as  t  — » -a.  s  — >  wo  and  w,  y,  z,  r  — >  o  (16) 

5  for  all  times  w  +  y  +  z  +  r  +  s  =  w  M  7) 

The  quantities  w.  y.  z.  r  and  s  from  equations  8  through  12  are  plotted  in  figure  2  versus 
time  for  a  particular  set  of  rate  constants.  These  are 

k,  =  1.0 
k.  =  0  5 
k,  =  0  25 
k4  -  0.125 

These  values  of  the  rate  constants  were  chosen  to  show  that  the  kinetic  model  predic¬ 
tions  can  be  similar  to  available  experimental  results.  The  curves  of  figure  1  are  actual 
measured  concentrations  of  2NDPA  and  its  derivatives  as  determined  by  thin  layer 
chromotography  as  a  function  of  time  at  90°  C  for  a  double  base  rocket  propellant  (ref 
8)  The  curves  of  figure  2  reproduce  qualitatively  many  of  the  features  of  figure  1 ,  such 
as  equations  13  through  17.  thus  justifying  assumptions  made  in  setting  up  the  kinetic 
model.  In  addition,  the  decay  of  the  concentration  of  2NDPA  in  a  double  base  rocket 
propellant  has  been  found  recently  to  be  exponential  with  a  finite  initial  slope  as  pre¬ 
dicted  by  equation  8,  the  initial  slopes  for  the  dinitro  DPAs  have  finite  values  as  pre¬ 
dicted  by  equation  9,  the  initial  slope  fortinitro  DPA  is  zero  as  predicted  by  equation  10, 
and  the  concentration  of  dinitro  DPA  has  a  maximum  as  a  function  of  time  as  predicted 

by  equation  9.12  The  isomers  for  the  model  were  choosen  on  the  bases  of  those 
reported  at  that  time  (ref  8).  More  recently,  other  isomers  have  been  reported,  mostly  at 

other  temperatures.1 2  It  is  straight  forward  to  include  all  isomers  in  the  model  in  such  a 
way  as  to  separate  the  concentrations  of  the  various  isomers. 

As  noted  above  it  is  assumed  in  writing  iiie  right  side  of  the  differential  equations  2 
through  6  that  the  concentration  of  diffusing  (NOx)  is  independent  of  time  for  the  times  of 

interest.  Thus,  the  rate  of  production  of  NOx  must  equal  the  rate  of  capture  of  NOx. 


1 .  G.  K.  Landram  and  R  J.  Bazuk,  Private  Communications 

2.  D  Robertson,  Private  Communication 
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This  is  at  best  only  an  approximation  and  may  account  for  some  of  the  differences 
between  figures  1  and  2.  The  model  does,  however,  indicate  how  diffusion  can  play  a 
role  in  the  thermal  degradation  of  propellant. 

DISCUSSION 

In  summary,  it  can  be  said  that  although  there  is  ample  evidence  for  diffusion  in 
the  thermal  degradation  of  NC  and  NC/NG  propellants,  there  is  little  in  the  way  of 
experiments  to  understand  this  diffusion.  This  thermal  degradation  may  be  rate  limited 
by  diffusion,  i.e.,  the  time  for  NOx  to  diffuse  from  the  site  of  production  to  the  site  of 

capture,  it  may  be  rate  limited  by  bond  breaking  to  produce  NOx,  it  may  be  rate  limited 
by  the  rate  of  reaction  of  NO  with  the  stabilizer  molecule  after  arrival  of  NO  at  the 

J  X  X 

stabilizer  site,  or  it  may  be  limited  by  some  combination  of  these  processes.  It  may  be 
possible  to  distinguish  between  some  of  these  possibilities  by  kinetic  studies  of  thermal 
degradation  of  NC  and  NG  with  and  without  stabilizer.  Global  activation  energies  may 
be  gotten  from  the  temperature  dependence  of  the  rate  constants  of  the  kinetic  model  fit 
to  experimental  data.  However,  additional  work  will  be  necessary  to  understand  the 
meaning  of  these  energies.  The  kinetic  model  does  give  a  qualitative  description  of  the 
conversion  of  stabilizer  into  its  derivatives  and  the  model  can  be  easily  extended  to 
include  all  of  the  derivatives,  as  noted  above.  It  should  be  possible  by  a  comparison  of 
the  model  with  experiment  to  determine  the  conditions  for  the  release  of  red  fumes. 
Thus,  the  model  could  then  be  used  to  predict  the  lifetimes  of  propellant  in  other  cases. 
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ALCULATIONS  FROM  KINETIC  MODEL 
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Figure  2.  Calculated  concentrations  of  2NDPA  and  its  nitro  derivatives  as  a  function  of  time 
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U  S.  Army  Armament  Munitions  and  Chemical  Command 
ATTN:  SMCCR-RSP-A 
Aberdeen  Proving  Ground,  MD  21010-5423 

Director 

Ballistic  Research  Laboratory 

ATTN:  AMXBR-OD-ST 

Aberdeen  Proving  Ground,  MD  21005-5066 

Chief 

Benet  Weapons  Laboratory,  CCAC 

Armament  Research,  Development  and  Engineering  Center 
U.S.  Army  Armament  Munitions  anc  Chemical  Command 
ATTN:  SMCAR-CCB-TL 
Watervliet,  NY  12189-5000 

Commander 

U.S.  Army  Armament,  Munitions  and  Chemical  Command 
ATTN:  SMCAR-ESP-L 
AMSMC-DS 

Rock  Island.  IL  61299-6000 
Director 

U.S.  Army  TRADOC  Systems  Analysis  Activity 
ATTN:  ATTA-SL 

White  Sand  Missile  Range,  NM  88002 
Director 

U.S.  Army  Material  Technology  Laboratory 
ATTN:  AMXMR-MMS 
AMXMR-H 

Watertown,  MA  02172 
Commander 

U.S.  Army  Materiel  Command 
ATTN.  AMCCN-C 

AMCLD,  D.  Vitali 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 
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Commander 

U  S.  Army  Armament,  Munitions  and  Chemical  Command 
ATTN:  AMSMC-DCG 
AMSMC-DSM-A 

Aberdeen  Proving  Ground,  MD  21010-5423 
Commander 

U.S.  Army  Depot  Systems  Command 
ATTN:  AMSDS-QV 
AMSDA-T 
AMSDS-SM-SPA 
AMSDS-CG 

Chambersburg,  PA  17201-4170 
Commander 

U.S.  Army  Toxic  and  Hazardous  Materials  Agency 
ATTN:  AMXTH-CO 
AMXTH-ES 
AMXTH-SE 
Edgewood  Area 

Aberdeen  Proving  Ground,  MD  21010-5423 
Commander 

U.S.  Amry  Western  Command 
ATTN:  APOP-NC 
APLG-MU 

Fort  Shatter,  HI  96858-5465 
Commander 

U.S.  Army  Central  Ammunition  Management  Office  (Pacific) 

ATTN:  SMCCA-CO 

Fort  Shatter,  HI  96858-5465 

Commander 

Johnston  Island  Chemical  Activity 
ATTN:  APGP-CP 
APGP-CM 
APGP-CH 
Johnston  Island 

APO  San  Francisco,  CA  96305 
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